Detection of noise-corrupted sinusoidal signals with Josephson junctions 
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Abstract 

We investigate the possibility of exploiting the speed and low noise features of Josephson junctions for 
detecting sinusoidal signals masked by Gaussian noise. We show that the escape time from the static 
locked state of a Josephson junction is very sensitive to a small periodic signal embedded in the noise, 
and therefore the analysis of the escape times can be employed to reveal the presence of the sinusoidal 
component. We propose and characterize two detection strategies: in the first the initial phase is supposedly 
unknown (incoherent strategy), while in the second the signal phase remains unknown but is fixed (coherent 
strategy). Our proposals are both suboptimal, with the linear filter being the optimal detection strategy, but 
they present some remarkable features, such as resonant activation, that make detection through Josephson 
junctions appealing in some special cases. 

PACS numbers: 07.05.Kf, 85.25.Cp, 84.40.Ua, 05.10.Gg 
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I. INTRODUCTION 



Detection of sinusoidal signals corrupted by Gaussian noise is, in principle, a completely solved 
problem. It can be proven that in the framework of statistical decision theory (Neyman-Pearson 
scheme) an adapted linear filter implemented via a Fourier transform is the optimal choice yj]. 
Unfortunately it is not always possible to apply such optimal detection strategy, due to a huge 
amount of data to process or to the extreme weakness of the signal to be detect. As an example 
we mention the all-sky all-frequency search for gravitational waves emitted by a pulsar [2], the 
search for continuous monochromatic signal in radio astronomy [3] and the detection of THz ra- 
diation ||4|]. When the optimal choice is not applicable, suboptimal strategies should be compared 
1 5]: for example it has been proposed to use generic bistable systems JaEzI- The general idea is 
the following: the oscillatory signal corrupted by noise is applied to a nonlinear bistable element 
that transforms the original signal to a new one, simpler to analyze. In other words, instead of 
applying a linear matched filter to the source signal, one inserts a nonlinear element; such inser- 
tion cannot improve the overall performance [8], but can make detection easier or can reduce the 
computational or memory burden. 

Among nonlinear elements (like generic bistable devices [6] and superconducting systems l^l) 
a good candidate is the Josephson Junction [10] (henceforth JJ) for two important reasons: 1) 
JJs are extremely fast elements, that can easily operate above 100GHz and up to the THz region 
llllll : 2) JJs , being superconductive elements, can be cooled close to absolute zero or to the 
quantum noise limit [120 , thus lessening the additional contribution of thermal noise from the 
bistable element. The detection of small signals in the presence of noise with JJs has been reported 
in two co nfig urations: investigating the magnetic flux trapped by a JJ closed in a superconducting 
loop JilLi (thus forming a SQUID llOJ) or analyzing the switching from the metastable locked 
solution to the running state of an underdamped isolated JJ [15]. The SQUID exhibits a standard 
double well potential equivalent to the classical case, well known in the context of stochastic 



180, although suboptimal 



resonance |16L Il7|] as a relevant phenomenon for signal detection 
[8J]. However SQUIDs, even if they are very sensitive and well studied systems, meet a principal 
drawback in their limited bandwidth, typically bounded to few KHz Ill8l . ll9ll . The escape from the 
static solution of single JJ II 1511 . being extremely fast, seems to be promising for signal processing. 
In this paper we aim to characterize such a device as a detector. 



The detection scheme that we propose is closely related to the general problem of passage over 
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a time-dependent barrier. More precisely, when the JJ passes the maximum of the time-dependent 
potential, switching from the static to the running solution, it sets the first passage time across an 
absorbing barrier [20Q. So the mean time between the initial state and the switch (or absorption) 
can be interpreted as a mean first passage time. Signal detection is then based on the possibility 
of ascertaining if the first passage times are due to a constant barrier (no signal present) or by an 
oscillating barrier (signal is present). 



II. MODEL AND DETECTION STRATEGIES 



The tunnel of Cooper pairs and normal electrons through a JJ biased with a sinusoidal signal 
of amplitude S , a noise signal S(t) and an additive thermal noise N(t) can be modeled by the 
Langevin equation lioll : 

+ WE + /cSin( ^ ) = Ib + SosinCm + <Po) + S(t) + N(t). (1) 
Eq. (0Q) is called underdamped because of the presence of the second derivative, or non negligible 
capacitance. In normalized units the equation reads: 

Cp + aip + sin(<p) = 7 + e sm(ut + (p ) + £(f) + 77(f). (2) 

Here a = (h/2eI c C) 1 ' 2 /R is the dissipation (I c , R and C are the critical current, the resistance 
and the capacitance of the junction), 7 = Ib/I c > tne dc bias due to the physical bias current 1b- 
These parameters can be adjusted in the experiments, although it is not easy to control dissipation 
(determined by the normal resistance R) one could in principle insert shunt resistors to achieve 
the desired level of dissipation. We have set a = 0.2 and 7 = 0.6 without performing any pa- 
rameter optimization to devise the best compromise between signal detection performances and 
experimental simplicity: we have just fixed the parameters for numerical convenience. Further- 
more e = S Q /I C is the amplitude of the ac signal, uij = [2el c / '(hC)} 1 ^ 2 the Josephson frequency, 
and t = uijt the normalized time. Overdot denote the derivative respect to t. 

In Eq. (0) two random terms appear, 77 that represents thermal current fluctuations of intensity < 
77(f) 77(F) >= 2a96(t — t'), (6 = 2eksT / (%I C ) is the normalized temperature) and£ that represents 
an additive noise with autocorrelation function of intensity D, viz < £(f)£(f') >= 2DS(t — t'), 
corrupting the external signal. To assume that the signal is only corrupted by an additive term 
is a simplification: noise can affect the signal in several ways, for instance as a bandpass noise, 
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multiplicative noise, phase noise, or as frequency fluctuations. We have focused on the simple 
case of additive noise (that is a paradigm in signal processing); however we expect that the results 
can be indicative of the behavior also for other noise sources. For instance frequency fluctuations 
of the drive in a JJ can be treated, in some limits, as an additive noise 112 ill . Quantum fluctuation 
contribute as an equivalent thermal source of temperature 9* = euj/(7rl c ) [|14l 12211. that can be 
confused with the stochastic effects [23]. Thus it is in principle always possible to decrease the 
temperature, proportional to 9, to reduce fluctuations to an unavoidable quantum noise level. 



A washboard potential is associated with Eq. © [1 



24] 



U(ip) = — •yip — cos(t^). 



The washboard potential © for 7 < 1 gives rise to a barrier [24] 



AU(j) = 2[y/l -7 2 -7 cos -1 (7)]. 



(3) 



(4) 



If the oscillating current is zero (e = 0) for low noise (9 « D, D « AU) escape occurs at a 



rate r 



r oc r K exp( 



■ AU, 
D ' 



(5) 



Such a rate is related to the average escape time r = l/r (rv is the Kramer prefactor, see Q20Q). 
The escape times can be directly measured in experiments fl23|,|25|]. In fact when an underdamped 
junction leaves the metastable zero voltage state it switches to a running state with which is asso- 
ciated a voltage. It is therefore possible to measure the time elapsed from the application of the 
signal to the escape. 

The possibility of experimentally determining the probability distribution of the escape times 
and the average escape rate is the key feature that we want to exploit for signal detection. More 
precisely, we propose to take advantage of the exponential dependence of the mean escape time 
upon the barrier height. The escape waiting time is highly sensitive to the amplitude of the signal 
feql as noticed since the pioneering experiments Il27ll that reported the striking changes into the 
escape time distributions for the resonant activation of JJ. For the case of no signal, the expected 
residence time distribution is exponential, while for oscillating barriers the distribution is modi- 
fied and the average changed. The phenomenon has been thoroughly analyzed for overdamped 
systems [28] with colored noise J30I1 and measured for underdamped systems [15]. For square 
wave signals, where a sudden switch of the barrier © between two values occurs, a remarkably 
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accurate estimate has been proposed Il26ll . Unfortunately, when a sinusoidal signal is embedded 



in noise and the barrier © is modulated, the formula of Ref. I126ll is no longer valid. Moreover, 



the method used in Refs. [[15 



2611 assumes that the signal is applied with a known phase (ip in 



Eq. © ) and therefore the escape time distribution is a function of the initial phase Ill5ll . In the 
context of signal detection, the frequency of the signal might be unknown, and therefore such ap- 
proach could be inapplicable. In this paper we propose to use two procedures. The first can be 
implemented without knowledge of the phase (incoherent detection), and only relies on detection 
of the escape time and subsequent reset of the system to the bottom of the potential well (see Fig. 
1). The second strategy (coherent detection) uses the phase parameter (that is supposed constant, 



although unknown) and hence is more similar to the approaches 115 



26J. 



The procedure that we propose for the incoherent detection strategy is the following: a) the 
signal-noise is applied with an unknown initial phase ip Q in Eq. ([2]); b) when an escape occurs in 
the evolution of Eq. © (the junctions reaches the metastable state U max ) the time Tj necessary 
for such an escape is recorded and the signal-noise is arrested or disconnected from the system; 
c) the JJ is reset to the bottom of the potential well U m i n ; d) the signal-noise is reapplied, i.e. Eq. 
(O is integrated adding ojt.i to the previous initial phase. We have numerically checked that after 
few iterations the memory of the initial phase ip Q is lost and that the long term distribution of the 
escape times is independent of the choice of the initial phase. It is therefore also possible to apply 
the incoherent strategy if the initial phase is unknown. 

In the second case (coherent detection strategy) one assumes that the frequency to of the signal 
is fixed. Under such an hypothesis, we adopt the following procedure (for square signals it corre- 
sponds to the procedure of Ref. [15]): a) the signal-noise is applied with an unknown initial phase 
ip to a J J that is at rest in the bottom of the potential well ©. The phase between the signal and 
the JJ is therefore fixed, but concealed; b) when the escape occurs, the escape time Tj to reach the 
metastable state U max is recorded and the signal-noise is only arrested after a time that is a multiple 
of the period 2n/co of the external radiation to guarantee that the same (unknown) initial phase ipo 
of the signal is retrieved; c) the JJ is reset to the bottom of the potential well U m i n to reproduce the 
same initial conditions as at point a); d) the signal is restarted, i.e. Eq. © is reapplied, and it will 
have with respect to the JJ the same initial phase (p . We remark that the phase difference between 
the JJ and the signal, ip , is frozen but unknown. To interpret the data it is necessary to reproduce 
the results for all values of ipo. In fact the distribution of the escape times depends upon the initial 
phase ip , as will be discussed in Sect. IV. 



5 



We remark on the main difference between the two strategies: when the system switches from 
the locked to the running state in the incoherent strategy we immediately stop the application of 
the noise-signal, while in the coherent strategy we let the signal run to retrieve the initial phase. 

For both strategies if the signal is digitally recorded, as it is the case for the all-sky-all-frequency 
search of gravitational waves emitted by a pulsar, a preliminary digital/analog conversion could be 
required to obtain a signal in the JJ frequency range. In this case it is possible to apply the recorded 
signal with a much faster time, in the GHz range or as fast as the electronics allows (while the real 
time signal might be slower) thus achieving a considerable speed up. 

Typical Complementary Cumulative Distribution Function (CCDF) of the escape times in pres- 
ence of a sinusoidal signal (e = 0.2, broken lines) and without signal (e = 0, solid line) have been 
simulated by numerical integration of Eq. © with the Euler method I129ll and are shown in Fig. 
|2] Such a Figure can also be interpreted as the CCDF for the passage times of a particle in a well 
over an oscillating barrier (broken lines) or over a constant barrier (solid line). 

From the data it is clear that the CCDF for the two strategies, coherent and incoherent, are very 
similar, with just a change in the slope and therefore a change in the average escape time. 



III. RESULTS FOR INCOHERENT DETECTION 

In Fig. [3]we show the estimated average escape times through U max (see Fig. [T]) vs to of a phase 
particle subject to an external drive and detected with the incoherent strategy (dot-dashed line of 
Fig. [2]). Around the normalized Josephson resonant frequency u = (1 — 7 2 ) 1 ' 4 the escape time 
is very sensitive to the external signal, in fact it exhibits a large dip, i.e. a pronounced deviation 
of the average escape time respect to the unperturbed value even for small s / \fD. (We recall 
that e/y/D is proportional to the signal to noise ratio, SNR). At low frequencies the deviations 
are less pronounced but still relevant, while for high frequencies one recovers the Kramer escape 
time of the unperturbed system [20]. It is clear that there is a wide range of frequencies (0, ojq) 
where the analysis of the estimated average escape times < r > s can give a clue to the presence 
of an external drive corrupted by noise, as they are different from the average of the escape times 
without the signal, r . 

To make such analysis more quantitative, in Fig. |4]we have adopted the Kumar-Carroll (K-C) 
index dxc fl] : 
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> I (r)s ~ (t)n I ,~ 

d-KC = , =, (o) 

y/i (^(r) s + o*{t) n ) 

where (t)s,{t)n are the estimated average escape time over a prescribed interval T, with and 
without signal, respectively. We have also denoted with o~(t)s,o~(t) n the corresponding estimated 
standard deviations. The K-C index © is related to the receiving operating characteristics (ROC) 
of the detector, and as a rule of thumb is well approximated by the ROC of a matched filter with 
signal to noise ratio equal to (Lkc PlQ . In Fig. @]it is evident, as expected, that dxc grows by 
increasing the observation time T. The statistical analysis of Fig. |4] confirms the result of Fig. [31 
i.e. the existence of a peak at the geometric resonance lu {dxc — 25 at T ~ 10 4 (27r/u;o)), and 
of an interesting region for lower frequencies (dxc — 5), while the method seems inapplicable 
at frequencies above coo (dxc — !)■ From Figs. [3] and |4] it is also clear that there is no evidence 
of stochastic resonance [17] due to resonant activation, or prominent signal detection at matching 
between the external drive frequency and the noise induced unperturbed escape rate, as reported for 
instance in Ref . [1 1 511 . In contrast to stochastic resonance the results of Figs. [3]and|4]only display a 
noise independent resonance at tu . This difference can be ascribed to the peculiar manner in which 
the external signal is applied. In the incoherent strategy the system looses memory of the phase 
of the signal at the passage through the absorbing barrier (see Sect. II). Instead in Ref. || 151] the 
system is reset after each switch - such reset corresponds to the coherent detection to be analyzed 
in Sect. IV. 

In Fig. [5]we show the K-C index for several values of the ratio e/y/D. The purpose is to 
emphasize the behavior for small SNR, when reliable stochastic simulations are prohibitive. The 
data reduction demonstrates that the K-C index decreases roughly as an inverse power law of the 
SNR; the best fit procedure gives an upper bound of ~ 1.5 for the exponent. Extrapolating the 
results, it is thus possible to infer the behavior for SNR's lower than those reported in the figure. 
We recall that these performances refer to a specific signal duration T = 2 10 5 (2tt /ujq) in FigJ5] 
The K-C index increases extending the detection time T, insofar as the K-C index is roughly 
proportional to a/T, see also Fig. @] Thus combining the numerical power estimated by Fig. |5] 
and the square root dependence of the K-C index captured by the data of Fig. |4] it is possible 
to predict that to keep a fixed dxc while lowering e / VD, the detection interval should increase 
as T ps (e/y/T))~ 3 . Consequently, by extending the detection time T one can also achieve the 
desired level of the K-C index for low values of SNR. As expected, the matched filter requires 
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shorter detection time ( T « (e/yD) ). 

IV. RESULTS FOR COHERENT DETECTION 

In Section II we discussed the possibility of detecting a signal at a predetermined frequency 
u> employing a coherent detection strategy. In such strategy the signal is always applied with the 
same initial phase. The analysis of the estimated average escape time in presence of a signal with 
initial phase cp = is reported in Fig. |6] The data show a remarkable dip at low frequency 
(to ~ 0.01), that is not present in the incoherent detection strategy, see Fig. |3j To compare the 
performances we analyze the K-C index © in Fig. |7J This analysis, as does the simpler analysis 
of Fig. [61 confirms the existence of an interesting region (d KC — 10) at low frequency (to ~ 0.01) 
that is not present in the incoherent strategy (see FigH]). 

In Fig. [8] we show the dependence of the K-C index as a function of the signal scaled amplitude 
e/y^D (compare to Fig. [5]for the incoherent analysis). This analysis leads to the same conclusion 
as for the incoherent strategy, i.e. by increasing the detection time T ps (e/\/D)~ 3 one can 
achieve the desired level of the K-C index. 

A second resonant condition for the coherent detection strategy occurs at the angular frequency: 

u = C(<p )(—), (7) 



TO 

where t is the escape time of Eq. © for e = 0. The factor C(cpo) reads 1/4 for ip = and is 
different from the frequency condition of the standard stochastic resonance [7]. This difference has 
been also reported in the experimental findings of Ref. [15, 2(]]. In Fig. |9]the frequency relation © 
is further elucidated: we show for different values of D the estimated escape time. From the figure 
it is clear that the largest deviation of the signal induced average escape time < r >s occurs at the 
matching condition described by Eq. © over a wide range of noise intensity, 0.04 < D < 0.2. 
Fig. |9]thus indicates the presence of a stochastic resonance between the noise (which determines 
To) and the frequency u of the external signal. We ascribe the change in Eq. © with respect to 
the traditional stochastic resonance to the peculiar choice of initial conditions. The reset of the 
initial conditions changes the effective waveform of the signal, see Fig. [Kjl Let us define P(ifo) 
as the time from the initial phase and the phase that corresponds to the maximum of the signal and 
therefore to the minimum of the barrier ©. Assuming that the effective waveform of the signal 
is of the type depicted in Fig. [TOl the matching condition for stochastic resonance shown in Fig. 

8 



m~l is that the barrier reaches a minimum in a time P(ipo) close to the noise induced escape from 
the lowest energy well. Such resonant condition depends upon the initial phase cp ; in Fig. [12] 
we report the scaled deviations of the average estimated escape time in presence of the signal as 
a function of the ratio 2ty/t lu for different values of the initial phase cp . We hypothesize that the 
horizontal axis positions of the relative minimal C(<^ ) -1 in Fig. [T2~a correspond to the resonant 
condition of FigQTJ In other words, we observe that a resonant condition occurs when the time 
tq(AU-) to escape the minimum barrier (denoted by ALL in FigfTTI) matches the time in which 
the oscillating barrier reaches a minimum: 



7o(AC/_) = P(<f ). (8) 

The resonance phenomenon described by Eq. ([8]) is evident in FigJT2a where the barrier reaches 
a minimum in a single ramp-up. Eq. ([8]) has been compared with numerical data in Fig. [13] and 
shows excellent agreement. This simple picture does not hold when the barrier non-monotonically 
reaches the minimum; in fact in FigJT2b the condition ([8]) is not valid for ip = it. 
We wish to underline the following remarks about Eq. ®: 

1) The prediction of Eq. ([8]) in the case cpo = 0, P((po) = (l/4)(27r/a;o), correctly describes 
the experimental finding of an harmonically driven underdamped JJ , see Ref . j26\ ; 



2) When the barrier decreases monotonically (— n/2 < (p < it/ 2) and the escape occurs in a 
single r amp-up of the signal (see Fig. [TOl) . the physical situation is similar to a fluctuating 



barrier [|26l 



32l-[34ll: the time in which the washboard potential reaches the minimum, P(<fo) 



in our notation, plays the role of the barrier flipping rate in Refs. 



126, 



32 



34] 



3) By varying the bias current 7 (which is an external parameter) one can set the most appro- 
priate barrier (see Eq.s ([4]) and © ) to match the resonant condition ©. 

A qualitative argument to understand the behavior when the signal does not reach a maximum 
in a single ramp-up (see FigfT2b) can be sketched following Ref. 
that a resonant activation occurs in two cases: 



3311 . where it has been noticed 



a) The low and high frequency limits of the escape times are identical: an extremum occurs if 
the derivative with respect to u of the barrier is nonzero; 

b) The low and high frequency limits of the escape times are different: an extremum occurs if 
the derivative with respect to to of the escape time is opposite to the difference of the limits. 



+ —7 r ■ (ID 



The low frequency limit of the escape time can be estimated with an heuristic argument (see 
the Appendix for details ): 

r(<p , u) « r K e b (1 + Awe cos(^ )) + O((ur(ip , u)) ), (9) 
where tk is the Kramer prefactor and the positive constant A reads: 

Tic AU(esin(u> n )+7) 

A = _Jl e J n) W ( £ sin(y,o)+ 7 ). (10) 

The high frequency limit (see the Appendix) is given by the time independent potential: 

r ac/( 7 ) 

T{(po, to) « r x exp ^— 

A consequence of Eqs. © and (fTTT) is that for < cp < n the low frequency limit of the barrier 
is lower than the high frequency limit, and therefore the escape time will be longer ( the reverse 
situation will be observed in the range — n < cp < ). From Eq. © we are able to compute the 
following derivative 

G?T(<A),W) Aj7( E Bin(y )+ 7 ) 

Z ~ T ^ 6 5 ^ecos(y>o) + 0(wr(yo,w)). (12) 

Eq. (fT2l) . valid in the low u; regime, shows that the corrections to the zero frequency limit 
r((po, 0) ~ r^e d have the same sign of the derivative of the external drive, i.e. 
sign(cos ipo). The results of FigfT2l confirm this picture: a resonance, or a non monotonic be- 
havior is observed for all initial phases but v?o = ±7r/2 (where the signal derivative vanishes and 
one should consider higher order corrections) and is most evident for ip Q = 0, n - where the signal 
derivative is at a maximum -. In the most effective cases , Lp = and cp = n, the former leads 
to the condition © with C((po) correctly predicted by Eq. © (see also FigfOT). while the latter 
((fo = 7r) exhibits a different frequency relation when the signal maximum not reached in a single 
ramp-up. 



V. CONCLUSIONS 



We employed the escape times of a JJ from the locked state as a statistic tool to detect sinusoidal 
signals corrupted by noise. We have determined the main features of such a statistical detector for 
both incoherent and coherent detection. For the former we found that the detector: i) is most 
sensitive to the signal in proximity to the junction resonance, ii) shows a clearly asymmetric 
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behavior in frequency: very low sensitivity above the resonance, and relatively good performance 
for lower frequencies. In the coherent case, which requires a priori knowledge of the applied 
frequency, the performance depends upon the initial phase. For an appropriate choice of the initial 
phase it is also evident a resonant activation, i.e. a matching condition between the noise induced 
escape time and the external frequency, Eq. ®, well fits the numerical data. Our analysis of the 
role of the initial phase extends the previously observed resonant escape Il26ll . Future research 
could be directed towards: i) development of better detection strategies; it) coupling two or more 
elements to exploit the properties of arrays of JJs [35] ; and Hi) a more accurate analytical treatment 
of phase dependent resonant activation phenomenon. 

We remark that the detection scheme based on JJs is suboptimal, but could prove fast and 
capable of operation at very low temperatures, with low intrinsic noise. It is therefore in niche 
applications where speed and reduced noise are essential that this approach could be considered. 
Finally, a few words of caution: while we have analyzed the physical principles underlying JJs as 
possible detectors, practical applications require a deeper analysis of circuit design and of techno- 
logical limitations. 

We wish to thank S. Pagano and I. M. Pinto for fruitful discussions and R. Newrock for a critical 
reading of the manuscript. 
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Appendix 



In this Appendix we give the asymptotic behavior of r(ip , uj) in the slow and fast frequency 
limit (u — > and u — > oo ) . To compute the slow frequency behavior of the average escape 
time we use, following Ref. |17l| . periodically modulated escape rates of Arrhenius type. In this 
connection we have: 

where 7(f) = 7 + 5 sm(ut + y? ) is the time dependent bias and t k is the Kramer escape rate pref- 
actor. Since we are interested in computing the asymptotic limit to — >■ 0, there exists a frequency 
such that tut <C w4(r) <C 1, where, due the exponential like tail of the escape time distribution, 
the probability of occurrence of an escape time greater than 4(r) is negligible. Accordingly, we 
expand the average escape time in Taylor series in the variable tut and due to the exponential like 
distribution we approximation have the estimates (t) « r((po,u)), (t 2 ) ~ 2r(<^ , w) 2 . The Taylor 
expansion of Eq. (Al) reads: 

r(v9 , u) « r K e { d <¥>0)+7) ^1 - ^r(v9 , w) cos(<p )A[/'(£ shi(y? ) + 7)) + O((a;r(v9o, u)) 2 ). 

(A2) 

In the previous equation r((p , tu) appear on the right side of the formula (A2), it can be consis- 
tently eliminated by an iterative substitution (and truncation) of (A2) in itself. By defining the 
parameter : 

A = --Ae » 0) 7) A£/ / ( £ sin(y ) + 7 ) (A3) 

Eq. (A2) can therefore be written: 

r(y? , w) « r K e b (1 + Awe cos(^ )) + O((wr(^ , ^)) 2 ). (A4) 

By using ©, the condition A > holds. 

The asymptotic limit u; — > 00 is obtained by noting that fast oscillation can not be followed by 
the system dynamic [20]. We have: 

AU(j) 



T(ip 0} u) « r^exp 



D 



+ 0[— . {Ah) 
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FIG. 1 : A picture of the potential barrier and of the detection procedure. When the phase (depicted as a 
dashed circle in figure) reaches the top of the barrier from the left, the system is restarted with a suitable 
initial phase from the bottom (black disk in figure) of the potential well U((p). 
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FIG. 2: The Complementary Cumulative Distribution Function (CCDF) of the escape time with signal 
e = 0.2 in the incoherent (dot-dashed line) or coherent ( dashed line) case. We also show the response 
to pure noise (solid line). The initial phase for the coherent strategy is ipo = 0. Parameters are: 7 = 0.5, 
a = 0.2, u = 0.86, D = 0.05, u = (1 - 7 2 ) 1/4 ~ 0.93 
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FIG. 3: Typical estimated (t)s as a function of u (T = 2 10 5 2tt/uo) for the incoherent detection strategy. 
Parameters are: 7 = 0.5, a = 0.2, D = 0.05, e = 0.1. The average escape time in absence of the signal is 
t = 134.6. 




0.001 



FIG. 4: The dependence of the K-C index as a function of angular frequency for parameters 7 = 0.5, 
a = 0.2, D = 0.05, e = 0.1 in the case of the incoherent detection strategy. The curves displayed are for 
different observation time T = 2ttN/ujq, where N = 10 p . 
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FIG. 5: The dependence of K-C index as a function of scaled signal amplitude e/y/D (proportional to 
the SNR) for the incoherent detection strategy.The parameters used are 7 = 0.5, a = 0.2, oj = ojq, 
T = 2 10 5 27T/W 



130 
120 
110 

A 

v 100 
90 
80 
70 

0.001 0.01 0.1 1 10 

a> 

FIG. 6: Typical estimated (t)s as a function of to (T = 2 10 5 2tt/ujq) for the coherent detection strategy. 
Parameters are 7 = 0.5, a = 0.2, D = 0.05, e = 0.1, cpo = 0. The average escape time in absence of the 
signal is tq = 134.6. 
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FIG. 7: The dependence of K-C index as a function of scaled angular velocity for the coherent detection 
strategy. Parameters are 7 = 0.5, a = 0.2, D = 0.05, e = 0.1, cpo = 0. The curves displayed are for 
different observation time T = 2-kN/ujq, where N = 10 p . 
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FIG. 8: The dependence of K-C index as a function of scaled signal amplitude e/y/D (proportional to 
the SNR) for the coherent detection strategy. The parameters used are 7 = 0.5, a = 0.2, u = ojq, 
T = 2 10 5 2vr/u; . 
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FIG. 9: Scaled estimated average escape time (t)s in presence of signal, vs scaled period 2tt/tquj for 
different values of D in the case of coherent strategy (initial phase <po = 0). The parameters used are 
7 = 0.5, a = 0.2, e = 0.1, tq is the e = escape time. The vertical grid line indicates the resonant 
condition ©. 
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FIG. 10: Effective waveforms for a system prepared with initial phase tpo = —tt/2,—tt/4, 0, it/A. The 
figures show the sinusoidal drive between the initial phase t^o and the maximum of the signal, tp = tt/2. 
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FIG. 11: The figures show the JJ potential U(tp) = — — cos(tp) as a function of the phase. We consider 
7 = 7 + ecos((^o) with initial phase ip$ = (top) or ipo — tt/2 (bottom); for illustrative purpose 7 = 0.5 
and e = 0.1. The time P(fo) is defined as the time between the application of the signal (with initial phase 
tpo = 0, — vr/2) and the lowest height of the barrier, AU-. The first passage time r is defined as the time 
to overcome the maximum of the barrier (the vertical dashed line). The resonance defined by Eq. © states 
that a minimum of the passage time occurs when are equal the time in which the external signal reduces the 
barrier to the lowest value and the escape time from such barrier, assumed to be static. 

20 



a) 



b) 




1 2 ~" 5 10 20 50 100 1 2 5 10 20 " 50 100 

2n/(T u>) 27t/(t w) 



FIG. 12: The dependence of the average escape time (t)s (scaled to To ) for the coherent detection strategy 
as a function of signal scaled period 2-k/uj for different values of the initial phase ipo . In panel a) we show 
the escape times when the effective waveform is monotonic, i.e. —tt/2 < (fo < tt/2; in panel b) we show 
the four limit cases 990 = 0, tt/2, tt, —tt/2. We underline that for ipo = ir the non monotonic behavior of 
the effective waveform ( the maximum is not reach in a single ramp-up) leads to a maximum in the escape 
time. The parameters used are 7 = 0.5, a = 0.2, e = 0.1, D = 0.07. 
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FIG. 13: The dependence of the resonant frequency as a function of 2n/P((po), proportional to the inverse 
time between the initial phase and the maximum of the signal, see Fig. QT] The dashed line corresponds to 
the prediction of Eq. ©. Parameters are 7 = 0.5, a = 0.2, e = 0.1, D = 0.07. 
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